A dynamical point of view of Quantum 
Information: Wigner measures 



A. Baraviera, C. F. Lardizabal, A. O. Lopes, and M. Terra Cunha 



Abstract We analyze a known version of the discrete Wigner function and 
some connections with Quantum Iterated Funcion Systems. 

Dynamics, Games and Science II, DYNA 2008, Edit. M. Peixoto, A. Pinto 
and D Rand, Springer Verlag (2011) 



1 Discrete Weyl relations 



This section fohows parts of [3]. Consider the Hilbert space % = C . Let 
{\k)}^~Q be an orthonormal base. Fix q;„, a^^ £ [0, 1] and define the foUowing 
matrices ?7jv, Vn G Mn{'C): 

N-l N-1 

t/jv := eT^*"" ^ eT^*'^»(A:|, \/jv := e^'"- ^ - 1| (l) 

k=0 k=0 

together with the identification |j) = \j mod N). Such operators are unitary 
and we have 
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C/jvIO =e^'{""+')|;), ywlO =e^'""|Z + l) (2) 

Defining n := (rii,n2) G Z^, we have that Un and Vn satisfy the discrete 
Weyl relations 

U^W;^' = e'^'^^^-VJi^U^' (3) 

Also, inspired in the continuous case, we define the discrete Weyl opera- 
tors: 

W^jvH :=e-''^"i"^t/;jiF^^ (4) 

Such operators satisfy 

W^in) = W{-n) (5) 

and 

WN{n)WN{m) = e*^''("'™)M^Ar(7i + m) (6) 
where cr(n,m) := nim2 — n2mi. 

When normalized, the discrete Weyl operators form an orthonormal base 
for M7v(C). In fact, using ^ and (g]), we have 



1=0 

N-l 

: ^ g-i-^(nin2+2ni(a„+i)-2n2Qt,)^/|; _j_ 



i=0 

= <5„,,o ^ e-^(""+0 ^ NS„,o (7) 
This allows us to obtain 

tr(W*N{n)WN{m)) = iV(5„,„ (8) 
and therefore for all A £ AIn{C), 

A=j^ Y.*'-{^Nin)A)wN{n) (9) 



where := {n = (ni,ri2) : < < iV - 1}. 
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2 Introduction to the Wigner function 

This section follows parts of [10] . Given a quantum system, we are interested 
in obtaining another form of representing the wave function ^{x). Such object 
will be the Wigner function, which will depend on two variables, moment 
and position. In order to understand such functions, we need to study the 
structure of phase spaces. 

The Wigner function consists of a special way of describing density opera- 
tors. In principle, we could say that density operators are a more fundamental 
structure than its Wigner representation. For instance, the Wigner represen- 
tation is unable to describe the density operators associated to two-level sys- 
tems. However, due to its simplicity, we will see that an understanding of the 
Wigner distribution gives us insight on certain aspects of density operators. 

Definition 1 Given a wave function ^{x), the Wigner distribution func- 
tion is 

1 /-oo 

W{q,p) = W^{q,p):=—J e^^^"' {q - -\^) {^q + -)ds (10) 
where above we are using Dirac notation 

{<i~^m = nq-'l) (11) 

{<Ir\q+l)^<P*(q+t) (12) 
Define the change of coordinates 

x = q+^^, = q- ^ (13) 

and then we obtain 

1 f°° 

Wiq,p)^— e*P(^-^ )(a;'|tf')(>fb)ds (14) 

That is, the Wigner distribution is obtained by calculating the product 
^{x')^* {x) and then applying the Fourier transform on s ^ x — x' . Such 
distribution has the following properties: 



/oo 
W{q,p)dp = {q\^){1'\q) =W{<1)? 
-OO 



(15) 
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W{q,p)dq={p\9){9\p) = \^{p)f (16) 

) 

/ W{q,p)dpdq = 1 (17) 

-oo J —OG 

where tp' is the moment representation of the wave function 'F. 

The Wigner function is real, but can assume negative or positive values. 
In this sense, it is not a density, but it is a kind of joint distribution of the 
position and momentum distributions. 

Now, note that (|14p can be written as 



W{q.p) = I e^P(^-^ '>{x'{\F){F\)x)ds = — I eiP^^-^ '>{x'\p\x)ds 

(18) 



27rhJ_^ ' VI /\ ly / 27rh 

where 



x = <7 + ^, x' = q-^ (19) 

where we define the density operator associated to a pure state as 

p-.^mii^i (20) 

The general definition for p includes pure and mixed states: 

P = Ep*I"^^)('^^I (21) 

i 

where Pi > and ^^Pi — 1. Such equation describes p as an incoherent 
superposition of pure state density operators !!?■,;) (tf^l, where 'Fi is a wave 
function, but not necessarily an energy eigenstatc. On equation pip the pi 
denote the probabilities of finding the system on the state jS',;). 

Hence, besides the usual probabilistic interpretation for finding a parti- 
cle described by a certain wave function at some position, we also have a 
probability distribution that such a particle can be found in different states. 





3 Discrete Wigner function 

This section follows parts of [7] and [11]. In dimension 1, the continuous 
Wigner function is in 1-1 correspondence with a density matrix p and is 
defined by 
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W,{q,p) = W{q,p) -.^—J j^^P/\q--\p\q+-)dX (22) 

Such function is uniquely defined by the following properties: [7]. [11]: 

1. W{q,p) e M 

2. If Pi and p2 are two density states then 

tr(pip2) = 271^1 J Wiiq,p)W2iq,p)dqdp (23) 

3. (Projection property) The integral along a line on phase space, described 
by aiq + a2P = 03, is the probability density that the measurement of the 
observable aiQ + gives 03 as a result. 



Remark Note that the Wigner function is always associated to a den- 
sity matrix. It would be more appropriate to use the notation Wp instead 
of W. When there is no possibility of confusion we will denote W. The pro- 
jection property stated above means, in other words, that the projection of 
the Wigner function along any direction of the phase space is equal to the 
probability distribution of a certain observable aiq + a2p, associated to that 
direction. Two special cases of this property are well-known: 



(24) 



W{q,p)dq 

is the probability distribution for the moment, and 

W{q,p)dp (25) 

is the probability distribution for position. For more details on these proper- 
ties, see pT] . 



We can write W as the expected value of a Fano operator, so we have 

Wiq,p)=tripA{q,p)) (26) 

where A can be written as 
1 



(27rfi)2 
1 



exp 



^{P - p) + ijiQ - q) dXdX' (27) 



(2..)^/ mA')e..4--(A',-Ap) 



dXdX' 



(28) 
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where 
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D(X,X') -.^exp --(Xp-X'Q) 
I n 



Also we can write A as 



Ml^P) = -^DRD* 



(29) 



(30) 



where above we write D = D{q,p) and R is an operator acting on positive 
eigenstates such that R\x) = | — a;). 



The proof that W satisfies properties 1 to 3 stated above follows from 
simple phase space properties. The fact that W{q,p) e K is a consequence of 
the fact that A{q,p) is hcrmitian. As for property 2, we can show that 

tr[A{q,p)A{q' ,p')) = ±^^5{q - q')5{p - p') (31) 

As a consequence, it is possible to invert equation (j26p so we can write 

p = 2nhj W{q,p)A{q,p)dqdp (32) 

Property 2 follows from the formula above. As for property 3, note that by 
integrating A{q^p) along a line on phase space gives us a projection operator. 
Therefore 

j S{aiq + a2p- a3)A{q,p)dqdp ^\a3){a3\ (33) 

where jaa) is an eigenstate of the operator aiQ + with eigenvalue 03. 
Later we will describe the proof of this property for the discrete case. 



Now we are interested in defining the Wigner function in the discrete case. 
The first step is to define a discrete phase space. Consider a Hilbert space of 
dimension N and define a base 

B, = {\n),n = 0,...,N~l}, 

which will be seen as a discrete position base. Now we define a base of 
moments 

Bp^{\k),k^O,...,N -1} 

A natural way of introducing the moment base from the position base is via 
the discrete Fourier transform. Then we can obtain the states of Bp from 
the states in B^ in the following way: 
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\k) = ^ ^cxp[27rmfc/7V]|n) (34) 

Therefore, as in the continuous case, position and moment are related by the 
Fourier transform. 



Remark We can relate the dimension of the Hilbert space with the Planck 
constant in the following way. We are supposing that the phase space has 
a finite area, which we can suppose equal to 1. In this area we can have N 
orthogonal states. If each state fills an area equal to 2TTh, we have N = 1/2ttH. 
So N plays the role of the inverse of the Planck constant and the limit as N 
goes to infinity can be seen as the semiclassical limit [7]. 



Given position and moment bases, we can define their respective displace- 
ment operators. For discrete systems, we can define translation operators U 
and 1^, in a way which is similar to what we have in ([T]) and ([2]), section [1] 

;7"|7i) |n + m), J7™|fc) ~ cxp[-27rimfc/Ar]|fc) (35) 

where the vector sums arc mod N. In a similar way the operator 1^ is a shift 
on moment base, and it is diagonal on positions: 

V"'\k) := \k + m), := exp[27rimn/iV] |n) (36) 

Then it is possible to show that 

ypf/q = g^^vqugyP, (37) 

the discrete Weyl relations ([3]) , seen on section[T] Let us also define a reflection 
operator as := | — n). We have that 

UR^RU-\ VR^RV-^ (38) 

The reflection operator is related to the Fourier transform in the following 
way. Denote by Uft the discrete Fourier transform, that is the operator 
whose entries on base are 

{n'\UFT\n) =exp[2Trinn'/N] (39) 

Then we have 

R = (40) 
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In order to define the discrete Wigner funetion, we still have to define a 
translation operator T and a point operator A, corresponding to the Fano 
operator defined in the continuous case. This is what we will do next. Define 

f{q,p) ■= IJ'^VP eyi^[iT:qp/N] (41) 

Such operators satisfy 

f(A<7,Ap) =f\g,p) (42) 

Remark In K.^ we define the translation operator with position q and 
moment p as 

f{q,p)=e-TMP~pQ) (43) 



Instead of definitions ([35]) and (p6| given for U and V we could, in principle, 
define U and V as the exponential of two operators Q and P, defined as being 
diagonal in Bx and Bp. However, infinitesimal operators Q and P satisfying 
the canonical commutation relations (CCR) cannot be defined over a discrete 
Hilbert space [4], [11]. Because of that we will use the finite cyclic shifts, given 
by (133 and 



Remark Due to technicalities, the phase-space can be taken to be a x iV 
or a 2N x 2N grid [7|. Typically we will be interested in phase spaces with 
even dimension and we will use the 2N x 2N grid (for instance, if = 2 the 
phase space has 16 points). Our following definitions will follow this choice 
as well. 



Let a = {q,p) be a point of the discrete phase space, with q and p assuming 
values between and 2N — 1 . Define 



A{a) y f (A,A')exp 27rz^^^^— Ml = J_u,rv-p^^^p,/n 

^ ' (2Ny ^ ' ' ^ i 2N \ 2N 

^ ' A,A'=0 

(44) 

We can express the translation operator in terms of A{a) by inverting the 
above definition and then we obtain the Fourier transform of A: 

f{n,k)= i(g,p)cxp[-i — (np-fcg)] (45) 

Note that as we defined the point operators over a lattice of 2N x 2N points, 
we get a total of AN"^ operators. However, such set is not independent. That 
is, we can show that 
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A{q + agN,p + CTpN) = i(q,p)(-l)'^p9+'^5P+'^,<^pA^ (46) 

for (Tq,ap ~ 0, 1. So we have that N'^ operators define the remaining ones. 
Define 

Gn -.^ {a = {q,p) : < q,p < N - 1} 
And the set G2N will denote the entire lattice of order 2N. 

A relation between A and T is the following: 

A(a)A{a ) = T{a~a) (47) 

By taking the trace of the above equation we get 

tT{A{a)A{a')) = -^^Niq' ~ q)SN{p' - P) (48) 

where a and a' are in Gn and 

N-l 

^^('^) '-^nT. e-'"'^'"/^ (49) 

ri=0 

is the periodic Dirac delta function, which is equal to zero unless q = 
mod N. 



Definition 2 The discrete Wigner function is 

Wia) = Wpia) tr{A{a)p) (50) 

where a G G2N ■ 

These 4iV^ values are not independent because in a similar way to what we 
have for the operator A, we have 

W{q + aqN.p + apN) = W{q,p){-lY^''+''^P+''^''''^ (51) 

for (Tg, CTp = 0, 1. As the operators A{q) form a complete set, we can write the 
density operator as a linear combination of the A{a). So we can show that 

p = AN W{a)A{a) ^ N W{a)A{a) (52) 

Remark It is possible to show that the discrete Wigner function defined 
above satisfies properties 1 to 3 stated in the beginning of this section. Prop- 
erty 1 is a consequence of the fact that A{q,p) are hermitian operators. Prop- 
erty 2 follows from the completeness of the set A{a), which allows us to show 
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that 

tr{p^P2)^N ^ W^{a)W2{a) (53) 

The proof of the third property requires a brief analysis of the lattice Gm 
and we refer the reader to section 21 for details. 

Conclusions We have defined the Wigner function for systems over a 
Hilbert space of dimension N < oo. The Wigner functions is defined as the 
expected value of the operator A{a) defined over the phase space given by 
equation (|44|) . The definition is such that W{a) G M. and is such that we can 
calculate the inner product between states and gives the correct marginal 
distributions along any line over the phase space, which is the lattice G2N 
with AN'^ points. Also, the values of M^(a) on the sublattice On are enough 
to determine W in the entire space. 





4 Calculating Wigner functions 



In order to calculate the Wigner function of a quantum state, we will use 
([35)). (|36p and (l44l) sc we can write W in the following convenient form: 



Lemma 1 

Proof In the following calculations, recall that the inner product is linear on 
the second variable. Wc have that 

W{q,p) = tr{Ap) = ^exp[iTrpq/N]tr{UmV-Pp) 

N-l ^ N-1 



^exp[inpq/N] J2{^T'\U'^RV-Pp\n) = ^ exp[i7rpq/N] {U~''ri\RV~P p\n) 

i=0 1=0 
^ N-1 ^ N-1 

— exp[i7rpq/N] ^ {n-q\RV^P p\72) = — cxp[i7rpq/N] ^ {q-n\V~P p\n) 

1=0 4=0 

" ^^^"P^-^^P^/^^ ^^^^'i ^ n)\p\n) = 



4=0 
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— exp[iTTpq/N] ^ exp[~2TTip{q - n)/N]{q ~ n\p\j 



1=0 

Also, note that 

nTpq/N-2Tnp{q-n)/N='-^{q-2{q-n)) = '-^{2n-q) = ^{n-q/2) 



Hence, 



1 2Trip 
Wiq,p) ^^"^{q- n\p\n) cxp[-^(n - q/2)] 

n—Q 



□ 



We believe there is a misprint in [7] in the expression corresponding to the 
W{q,p) described by the claim of the above lemma. 

One can ask how Wp changes with p. Suppose first p is a projector from a 
wave ■0 which has norm 1 in C^. Suppose {a4>i+b(j>2) = V'l where -0, 02 have 
norm 1, and p — \ip >< ■01- Then, W^^ ^ aW^^ + bW^^. The linearity occurs 
only when p = X^i^il* that is, when p is diagonal. This in general 

do not happen for operators ji/; >< ■01 induced by a wave if). However, if 
p = {api + bp2), where p, pi,p2 are density matrices, then Wp = aWp^ + bWp^ . 

Example 1 Let N ~ 2, and let |?/') = a|0) + be a state superposition. 
Let Wi{a) and ^2(0?) be the Wigner functions for |0) and respectively. 
We have that the Wigner function W for lip) is such that 

Wia) = \a\^Wi{a) + \b\^W2{a) + 2i?e{a6*(l| A(a)|0)} (55) 

In fact, note that 

W{a) = tr{A{a)p) = tr(A{a){\a\''\0){0\ + + a6*|0)(l| + a*6|l)(0|) 



\a\'^Wi{a) + \b\'^W2{a) + ab*tr{A{a)\Q){l\) + a*btr{A{a)\l){Q\) 
|opWi(a) + \b\'^W2{a) + ab*tr{{l\A{a)\{))) + a*btr{{Q\A{a)\l)) 



so the result follows. 







Let us remark a few properties of the Wigner function for a pure state 
p. In this case by expanding p in terms of the phase space operators as in 
equation (j52p and by imposing the condition p^ ~ p, we get 

W{a)=AN'^ ^ r(a,/3,7)W'(/3)VF(7) (56) 

/3,76Gjv 
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where the function r{a,P,j), which depends on 3 points (i.e., a triangle) is 
given by 



1 



ria, /3, 7) tr(A(a)A(/3)A(7)) = — exp — 5(a, /3, 7) , (57) 



2ni 



N 



of 2 of the 3 point (a, /3, 7) contain even q and p coordinates. Otherwise we 
define 

r(a,/3,7) :=0, (58) 

and in the above expression, vahd for even N, the value S{a, /?, 7) is the area 
of the triangle formed by these points (measured in units of the elementary 
triangle formed by 3 points which are one position apart from each other). 



Now we calculate the Wigner function for a position eigenvalue 

P9o = ko)(9o| (59) 
We obtain the following closed expression for W: 

= ^S^iq - 2go)(-l)^t(«-''°^ "^"'^ ^1 (60) 

We can also write the Wigner function of a state which is a linear superpo- 
sition: 

^ ±{\q,) + e-^^\q,)) (61) 
Again, we can obtain a closed expression for W, which is 

W{q,p) = i {W,, {q,p) + W,, {q,p) + AW,,,,, (g,p)) (62) 
where the interference term is 

AW,,,,,, {q,p) := l6^{q){-iyr> cos {^p + ^) (63) 

where 

2N 

q = qo + qi - A = (64) 

90 - qi 

This is an explicit expression for the calculation seen in example [TJ 



Now we make a few considerations on the time evolution of c^uantum sys- 
tems on phase space. If U is the unitary operator which determines the evo- 
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lution of a state, then the associated density matrix evolves in the foUowing 
way, 

pit+1) = Up{t)U* (65) 

By this fact, we can show that the Wigner function evolves in the following 
way: 

Wia,t + 1)^ Z^pW{l3,t) (66) 

where the matrix Z^p is defined as 

Zap Ntr[A{a)UA{l3)U*^ (67) 

Therefore the time evolution in phase space is represented by a linear trans- 
formation, which is a consequence of Schrodinger's equation. The unitarity 
imposes a few restrictions on the matrix Zap. In fact, since purity of states is 
preserved, the time evolution has to preserve the restriction given by equa- 
tion (|56p . Therefore, the matrix has to leave the function _r(a, /3, 7) invariant, 
that is, 

r{a\ fi\ 7') = Za'aZp>i3Zy^r{a, /3, 7) (68) 

The real matrix Z^/s contains all the information on the time evolution of the 
system. In general, such matrix relates a point a with several other points /3. 
So the evolution will be, in general, nonlocal, a unique property of quantum 
mechanics. In classical systems, the value of the classical distribution function 
W{a,t + 1) is equal to the value W{/3,t) for some point /3, which consists of 
a well defined function of a and t. However, we have in [7| a few examples of 
unitary operators which generate a local dynamical evolution on the phase 
space. 



To conclude this section, we calculate the Wigner function for a quan- 
tum channel A, as the ones considered for our analysis of QIFS. This is a 
straightforward calculation. Let Vi be linear operators, i = 1, . . . , fc such that 
V*Vi = I. Then A{p) = V,pV* G Mn- Hence, 

1 ^"^ r27ri 
WAip){q,p) ^ — ^((?-n|vl(p)|n)exp —p{n-q/2) 



2N 

n=0 



1 ^"^ rS> ■ 

1^ E E(« - ^\V^PVnn) exp \^p{n - q/2) 



2N 

n=0 1=1 
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N-1 k 

U(n - n^V^n\V*(■n^^ exn 

N 



^ E E(('^ - ^)np\V:[n)) exp -^p{n - q/2) (69) 



2iV 

,7V-1 



1 ^ r2^7 



2N 

n.j—0 i—1 



Writing p = J2j=o Pib'XiL Ej Pi = 1, we get 

' nJ(n ~ nWAi)(i\V*(n')) ran 

N 

(70) 

Therefore the Wigner function of A{p) is obtained in a simple way from the 
function for p. 



5 some properties of the discrete Wigner function 

We have seen in section [3] that the discrete Wigner function 

W{a) = tr{A{a)p) (71) 

satisfies properties 1 and 2. Now let us prove property 3. Let p = EiPiK)(*l; 
J2i Pi = 1 be a density operator. Denote by 

= {|n),n = 0,...,iV-l}, 

a position basis and 

Bp^{\k),k^O,...,N-l} 
a moment basis, as before, where 

\k) = ^ ^exp[27rmfc/^]|7i) (72) 



N 

To prove property 3, wc must show that as we sum the operators A(q,p) over 
the point of the phase space which lie over a line L, we obtain a projection 
operator. This implies that by summing the values of the Wigner function 
over all the points of a line we get a positive number, which can be interpreted 
as a probability. 



We begin by defining a line on the phase space. A line L is a set of point 
of the lattice, defined as 
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L = L(ni, 712, ^13) = {{q,p) e G2N ■ nip - = ng, < n, < 27V - 1} (73) 

Also, we say that two lines as parallel if they are parameterized by the same 
integers ni and 712. 

Now, let us show that as we sum the point operators A over a line, we get 
projection operators. So we are interested in the operator 

(9,p)6-L 



Since 5N{q) ~ jj Y^n=o ^ STrign/Af^ rewrite such operator as 

= ^ A{q,p)52N{nip-n2q~n^) 



q,p=0 
2N~1 2Ar-l 



2N ^ ^ 27V 

A=0 q,p=0 

1 ^^"^ 27r 
= ^ E ^'('^i''^2)exp[i— TisA] (75) 

A=0 

where we use the Fourier transform of A to obtain the last equality. Since 
T is unitary, we have N eigenvectors \4>j) with eigenvalues exp[— 27ri0j/iV]. 
Besides, such operator is cyclic and satisfies = /. Therefore as its eigen- 
values are N ^ th roots of unity, the are integers. So we can rewrite (j75p 
as 

1 ^ 2n 

= ^ E E«^P[-*^(2(/., - 773)A]|0,)(0,| 
A=0 j=0 
N 

= Y,^2N{2cl)j~nsM,){^,\ (76) 

Hence we have that Al is a projection operator over a subspace generated 
by a subset of eigenvectors of the translation operator T(7ii, 712). 



Example 1 For a line Lq defined by q = 713 (that is, ni = 1, 772 = OJ, the 
Wigner function summed over all point of Lq is 



E Wp{q,p) =Y,Wp{n3,p) = {n3/2\p\n3/2) (77) 
{q,p)eL, p 
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if is even, and equal to zero otherwise. 

More precisely, we have the foUowing proposition: 
Proposition 1 Let N be even and let p be a density operator. Then 

2N-1 

J2 Wp{2q,p)^{q\p\q), g = 0, 2, . . . , iV - 1 



and 



p=0 



2N-1 



Wp{2q + l,p)^0, g-0,2,...,7V-l 



p=0 



Proof First, to sec why the case q odd implies that the Wigner function 
equals zero, consider the expression for W given by 

Write p = Cj IjXjI, c, > 0. Then 

{q-n\p\n) = ^ Cj(g - n|j)(j|?i) (79) 
j 

which is 7^ if and only li j = q — 7i ~ n for some j. In particular, in order 
to have a nonzero inner product above, we must have that q is even, because 
q — n = n implies q = In. 

Now suppose that q = 2qQ. By the analysis above, we see that in the sum 
of the terms forming the Wigner function (eq. ([75| ). we only have to sum the 
indices such that the equation 

q — n = n 4^ 2go — n ^ n (80) 

is satisfied (recall that all calculations are made modulo N). Such equation 
has two solutions, namely n = qo and n = qo + N/2. To see that there are no 
other solutions for (j80|) . we proceed in the following way. From 2qo — n ^ n 
we get 2{qo — n) = 0. We know that n = and n = qo + N/2 are solutions. 
Also, note that a; = and x = N/2 arc solutions of 2x = 0. Now, if ?/ is a 
solution of 2a:: = then y — N/2 also is. Clearly if y is an element between 
and N/2 then 2y will be at most equal to 2N — 2, hence 2y 7^ 0. Finally, let y 
be an element between N/2 and N and by contradiction suppose that 2y = 0. 
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Then by the remark above we have that z = 2y ~ N/2 is also a solution and 
z is between and N /2. But there are no solutions for 2a; = between and 
N/2. This shows that 2x = admits only the solutions stated above. 

Now note that if n equals go then 

cxp ^P{n - g/2)] - 1 (81) 

If n = go + N /2, we have that the exponential above is equal to ±1, being 
positive or negative if p is even or odd, respectively. Therefore, for N even 
and q = 2qo, we have 

Wpi2qo,p) = ^i{qo\p\qo) ± (go + N/2\p\qo + N/2)) (82) 

where the sign ± depends on p. For fixed q and considering all possible p 
(i.e., p = 0, . . . , 2N — 1), we have that the second inner product above will 
have a plus sign in front of it in the N possibilities in which p is even and 
will have a negative sign in the N remaining possibilities. So 

Y,Wp{2qo,p) ^ {qo\p\qo) (83) 
p 

This concludes the proof. 

□ 

Corollary 1 If q is odd then Wp{q,p) = 0, for any p and any p density 
operator. 



Proof Follows from the first paragraph of the proof above. 
Definition 3 Let ip be a state. The VF-transform of ip is 

2JV-1 



□ 



</>(p) ~ J2 W^il^'^P) (84) 

q=0 



forp^0,...,2N-l. 



Let (j> be the W-transform of i/', and let J-ip be the discrete Fourier trans- 
form of ijj. 



Question: 
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\{Fi;){p)\^ = P = 0,1,...,7V-1 (85) 

Answer For N = 2 and ip = |0) or |1), the answer is yes. In fact, let 
1^) = |0) = (1,0). Then 



And 



Also 



And 



(J^|0))(0) = 4=^l(-^|0))(0)|^ = i 



</)(0) = EM^|o)(g,0) = J + + l + = ^ 



(^|0))(l) = i=^|(J-|0))(l)p = i 



</)(!) = E M^|o)(<7, 2) = i + + i + = i 



9 

Therefore in this case 



|(J-^)(p)|2 = 0(p), p^O,l (86) 



Now let IV') = |1) = (0, 1). Then 
■^11) = ^E^^P[2^*-?'/2]|j) = ^(|0> +cxp [27rz/2]|l)) = i=(|0) - |1)) 



And 



Also 



And 



Therefore 



0(o) = E^ii>('?'O) = i + o + l + o = J 



imm = -^^\{mm\'-l 



V)(l)=Ew^li>('?,2) = ^+0+i 

9 

|(J-^)(p)|2 = 0(p), p = 0,l (87) 
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Now let us write an example in which the state considered is mixed. Let 
V' = 1/V2(|0) + |1)). Then 



i=(|0) + |l)) + -i=(|0)-|l))] =|0) 



Then |(J"V)(0)P = 1 e = 0. Now let us calculate 0(p), p = 0, 1. 

By definition, we have = J2q ^^'il^ 2p). We can use the expression ([5^ : 



'AQ^ 0) - i (VK|0) {q, 0) + W^iD (q, 0) + Aosiq, 0) 
-^{q, 2) = i {W\o) (g, 2) + M^iD (q, 2) + AoAq, 2) 



Then 



Wv,(0,0) = 


2^ 


1 


f 0) 


W0(1,O) = 




-04 


1' 


M/^(2,0) = 




1 


f 0) 


2^ 




W^V'(3,0) = 


1(0 H 


-04 





1 

4 
1 
4 
1 
4 
1 
4 



which implies 0(0) = 1 = |(J-'7/')(0)p. Similarly, 



W^^(0,2) = i(i + i + + 0) = i 

W^v>(1.2) = i(0 + 0-i) = -i 
VF^(2,2) = l(i + l + + 0) = l 
W^(3,2) = 1(0 + 0- i) = -i 



and so (^(1) = = |(-FV')(1) 



(89) 
(90) 







Inspired in the calculation above, we prove the following lemma, valid for 
pure states only. After that, we will prove the result for density operators. 



Lemma 2 Let i/' = \m) e {|0), ...,\N~ 1)}, N even. Then 
|(J-V')(p)|' = 0(p), p = 0,l,...,iV-l 



(91) 
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Proof We have 

Af-l 



■^'"^^ 7^ ^ cxp[27rijm/iV]|j) 

So 



(J-|m))(p) = exp [2T:ipm/N] ^ |(J-|m))(p)|2 = 1 



Let us calculate = S5=o ^ ^|m) (^7 2p)- By the corollary[Tl wc only have 
to sum the even q. Then = X)^!^^ 2p). By proposition [1] we 

get, using expression ([5^ . that 

^^^(2^0,^) = ^((golpko) + (go + N/2\p\q„ + N/2)) (92) 

where the sign of the second inner product is positive because 2p is even. 
Now note that only one of the inner products above can be nonzero, because 
p is pure, by assumption. Moreover, p pure implies that such inner products 
are equal to 1. Finally, since q varies between and 2N — 1 we have exactly 
two nonzero terms in the sum of 4'{p) namely, the terms corresponding to the 
m and m + N /2 indices. Hence, 

<t){p) = l/2iV + l/2iV = l/N = |(J-|to))(p)|2 

This concludes the proof. 

□ 

The following result, inspired in the previous one, completes proposition[Tl 
which related the discrete Wigner function with the base of position vectors. 
Now we do the corresponding work for the basis of momentum vectors. 

Proposition 2 Let N be even and let p he a density operator. Let \p) he 
a vector of the momentum hasis, that is, obtained via the discrete Fourier 
transform of a position base vector: 



Then 



2N-1 



J2 Wp{q,2p)^{p\p\p), p = 0,l,...iV-l (94) 

2N-1 

J2 Wp{q, 2p + 1) = 0, p = 0, 1, . . . iV - 1 (95) 



g=0 
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Proof Let us calculate (f){p) ~ Y^'^=o^ ^p{Qi 2p). By corollary[Tl we only have 
to sum the even q indices. Then (j){p) — X]^o^ M/p(2'7, 2p). By proposition [T] 
we get, using expression (j82|) . that 

Wpi2q, 2p) = ^{{q\p\q) + {q + N/2\p\q + N/2)) (96) 

where the sign of the second inner product is a plus because 2p is even. Write 
p = Ci|i)(z|. Take, for instance, q = 0. Then 

Wp{0, 2p) = ^((0|p|0) + (0 + N/2\p\0 + N/2)) 
= ^C£cMi)m + {N/2\i)m/2)) = ^(co + c^/2) (97) 

i 

As we know, Wp(l, 2p) = 0. Take g = 2, then 

W^p(2,2p) = ^(ci+cjv/2+i) (98) 

and so on (noting that wc always have zeroes when q is odd). In this way, we 
end up summing all Ci coefficients Ci twice (because q varies between and 
27V — 1) and we get that 

2N~1 ^ ^ 

cj>{p) = Wp{q, 2p) = -(co + ci + • • • + C2N^i) = ^ (99) 

9=0 

By the calculation above, we only have to calculate {p\p\p) and show that 
such number equals 1/N. Recall that the inner product we consider is linear 
on the second variable, so we write p = Cm\m){'m\ and then: 



^ N-l N-l 
{P\P\P) = cxp [-27ryp/iV] cxp[2T:ilp/N]{j\m){m\l) 



j=0 1=0 

""^^"'J^Y'^^P [-27rijp/7V] cxp [2mmp/N]{j\m) ^ ^™ = 

m j—0 m 

(100) 

□ 

Conclusion By propositions [T] and [2] we have for the discrete Wigncr trans- 
form that if N is even and p is a density operator then 
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2N-1 2Af-l 

Y,Wp{2q,p) = {q\p\q), ^ VF,(2g+l, p) = 0, = 0, 1, . . . , iV-1 (101) 

p=0 



p=0 

and if 



then 

2JV-1 



N-l 



^ ^cxp[27rijp/7V]|j) 



(102) 



3=0 



2N-1 



Wp{q,2p)^{p\p\p), J2 Wpiq,2p+l)=0, p = 0,l,...iV-l (103) 



q=0 



q=0 



Such expressions are the discrete analog of the result we have for the contin- 
uous Wigner function, namely the result that relates the marginals with the 
Fourier transform J-": if p = |V')(^| then 



Wp{q,p)dp=\i'iq)\^, 
See [S] for more details. 



Wp{q,p)dq^\Ti'{p)f' 



(104) 







Example 2 Denote by Wp the matrix with entries Wp{q,p) for q,p = 
0, ...,2A^-1. For instance, if N ^ 2 and writing \0) = (1,0) and \l) ^ (0,1), 
we have that Wp contains the image of the Wigner function for each point of 
the phase space. We immediately notice that the integral over all space equals 
1: 

.1111 



W|o)(o| = 



4 4 4 4 



4 4 4 4 

Vo / 



, W|i)(i| 



/- --- --\ 

I 4 4 4 4 \ 



1 1 1 1 

4 4 4 4 

Vo / 



(105) 







Example 3 Denote by Wp the matrix with entries Wp{q,p) for q,p ~ 
0,...,2N-l. Let N = 4,, and writing |0) = (1,0,0,0), |1) = (0,1,0,0), 
|2) = (0,0,1,0), |3) = (0,0,0,1), we have, in a similar way as seen in the 
previous example, that 
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W|o)(o| = 



|i>(i| 



W|2)(2 



1 


1 


1 


1 


1 


1 
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VV|3){3| 
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00000000 
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00000000 
00000000 
00000000 
iiiiiiii 

88888888 

VO 0600000/ 



Remark 1 What occurs in general for pure states: the Wigner function 
^l9o>(9ol except in two hnes, located in q= 2{modN). When q = 2qo, 

W assumes the value 1/2N, and when q ~ 2go i N, W assumes the value 
1/2N for even values of p and ~1/2N for odd values. Such oscillations are 
typical of interference fringes and can be interpreted as arising from the 
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interference between the line q ~ 2qQ and a mirror image formed at a distance 
of 2N from 2qQ, induced by the periodic boundary conditions [7]. 

Remark 2 The fact that the Wigner function assumes negative values 
in the interference line is essential for one to be able to recover the correct 
marginal distributions. Summing the values W{q,p) along a vertical line gives 
us the probability of measuring g/2, which should be equal to 1 if g = 2qo, 
and equal to zero, otherwise. 



A natural question is to try to understand the action of the operator which 
defines QIFS in the dual variables p. This is the purpose of the next results. 

Lemma 3 Let A(p) = VipV* and define F{p) ~ J-pJ-*, where J- is any 
unitary map. Then there is G : A4n ~^ Mn such that the above diagram 
commutes: 



Mn > Mn 



A 



F 

Mn > Mn 



(106) 



Proof First, note that F ^(p) = T* pT . Also T is unitary, therefore we have 
J"-i = J"*. Define G = F o Ao F^\ Explicitly, 



Gip) = FiJ2 v,^*p:fv*) = :f[J2v,^*p:fv: 

i i 

i i 

where Vi = TViT* . And a simple inspection shows that 
F{A{p))^G{F{p))=Y,TV,pV:T* 



□ 



Example 2 Consider N = 2. Then the discrete Fourier transform is given 
by 



In this case we have J- ^ = T . Let 
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where the pij form a column stochastic matrix P. Then lemma [5| for this 
example shows that G{p) — VipV* , where 



2"' \-l-lJ' 2" V-1 1 

Then, from pn + P2i ~ 1. P12 + P22 ~ 1 md writing 

_ / pii P12 

\P21 1 - Pll 

we get from f/ie expression 

f{A{p))^g{f{p))^Y.^v.p^:f* 



\ PllPll +Pi2(l - Pu) - 5 



PllPll +Pl2(l - Pll) - 5 ^ 



(110) 



2 2 

7n the case that the vector tt = (pii, 1 — pii) is fixed for the stochastic matrix 
P , we can rewrite the expression above as 

F(^(p))=G(F(p))= f 5 /"-l") (HI) 



Lemma 4 Define A : Mn ^ Mn, A{p) = J2i^iPW ^ w**'^ ^» linear, 
Si ^r^i ~ '^'^'^ '^^(p) associated discrete Wigner function. Then 

given {q,p) there are Mi = Mi{q,p) such that 



WAip){q,p)=Y.tr{AUpM*) 



Proof First, as A{q,p) is hermitian, we have a decomposition 

A = UDU^^ 



where U is unitary and D is diagonal (and real). Then 

^1/2 ^ 
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where {A^/'^Y = ^> D^/^ is the diagonal matrix whose entries are the positive 
square roots of the entries of D. Then 

WA(p){q,p) = tr{A{q,p)A{p)) = tr(i^K:P^;) = Y,tr{AV,pV*) 

i i 

= J2 triA^^''V,pV*A^^'^) = J2 tr{UD^/^U-^V,pV*UD^/^U~^) (112) 

i i 

Defining M, = UD^/^Lf-'^Vi and noting that JJ-^ = U*, we can write 

i 

□ 
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